On Peres approach to Pradkin-Bacry-Ruegg-Souriau's perihelion vector 
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We solve explicitely the differential system obtained by Peres for the construction of a conserved 
vector associated to any central potential. We then obtain a very direct access to the discontinuous 
behavior of this Fradkin-Bacry-Ruegg-Souriau perihelion vector. 
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I. INTRODUCTION 



Since the pioneer works of Ermanno, Bcrnouilli and Laplace during eighteenth centurji^i^i^ we know that the Kepler 
' system presents a specific vector conserved quantity, the so-called Laplace-Runge-Lenz vector, which provides in 
, particular a very simple access to the orbit equation. Today, the existence of generalized Laplace-Runge-Lenz vectors 

for various extensions of the Kepler problem have been put in evidence^. 
' This supplementary constant of motion has a dynamical origin which is directly linked to the existence, for these 
Qh' systems, of a symmetry group larger than the space-time's one^i^i^. As established by Jauch and Hill^ in the 2D 
case and by FradkiniS in the 3D case, the isotropic harmonic oscillator shares the same feature and admits also a 
supplementary dynamical conserved quantity but of tensorial type, the Fradkin-Jauch-Hill tensor^'^'^^. 
j3 ' Bacry, Ruegg, Souriau^iS and Fradkin^- made an important additional step when they showed that all three dimen- 
sional dynamical problems involving central potentials possess the extended symmetry algebra O4 and SU3, a result 
soon generalized by Mukunda^^ for any hamiltonian system with n degrees of freedom. As noted by Bacry, Ruegg 
and Souriau in their seminal paper—, the general ability to construct such a local Lie algebra of dynamical symmetry 
^ I is not surprising since all the 2n-dimensional symplectic manifolds are locally isomorphic^^. However the following 
' step to ensure a global dynamical symmetry is to verify that the finite canonical transformations generated by this 
00 1 algebra form a group, a feature of the Kepler problem. As established by Stehle and Hani^, the existence of such a 
C ' "global" higher symmetry group, admitting finitely multivalued realizations on the phase space of the system, is a 
, necessary condition for the existence of degeneracies. 

Bacry, Ruegg, Souriau"^^, as much as Fradkii*i^, gave an explicit construction scheme for a vector with constancy 
properties associated to every central potential, Fradkin's derivation being very direct and explicit. Bacry, Ruegg and 
QP ' Souriau already noted that such a vector is exceptionally one-valued, as in the Kepler case. The first explicit demon- 
stration of the multivalued behavior of these Fradkin-Bacry-Ruegg-Souriau (FBRS) vectors have been established by 
Serebrennikov, Shabad^, Buch and Denmarti^ when these authors pointed out that in specific cases Fradkin vectors 
are in fact only piecewise conserved. 

In 1979, Peresi^, using a distinct approach, rediscovered the FBRS generalized Laplace-Runge-Lenz vector. He 
' obtained a differential system for the coefficients of the vectors which presents singularities at the apsidal positions. 
Without solving explicitely the differential system, Peres infered from this singular behavior the existence of discontin- 
uous changes of direction of the vector which is then only piecewise constant. The correspondence between Fradkin's 
and Peres approaches have been studied by Yoshida^'^. 

A complete answer to the question concerning the status of conserved quantity which could be attributed to the 
FBRS vector has been given by Holas and Marcb^. They showed explicitely that at each pericenter (or apocenter, 
depending on the chosen initial conditions), the FBRS vector changes its direction abruptly from an angle double of 
the apsidal one. A last attempt to obtain a "true" vector constant of motion in every conservative central force field 
has been made by Yan32i but, as shown a short time after by Holas and March2^, Yan's construction coincides in fact 
with the FBRS perihelion vector and then presents the same discontinuities. The generalized Laplace-Runge-Lenz 
vector can also be deduce in a direct and elegant manner from the equation of motion-'^. 

Note finally that, in a very recent contribution, Ballesteros and al.— use FBRS type vectors to establish an optimal 
extension of Bertrand's theore m^^i^^ to curved spaces. 

In this paper, we give an explicit solution to the differential system obtained by Peres^^. Using then a complex 
formulation, particularly convenient for the treatment of planar problems^^i^, we recover in a very simple manner 
the discontinuous behavior of the FBRS vector studied by Holas and Marcb^. 
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II. BASICS OF 2D MOTION IN A CENTRAL POTENTIAL 

1. Complex formulation 

We consider a planar motion V [t) = (^[t]) ^ ^ for a particle of mass m = 1 submitted to a potential U{'V), 

eventually singular at the origin. We choose to adopt a complex formulation where we represent the position by its 
corresponding affix z[t) = x[t) +iy(t), the potential being then viewed as a real valued function of z, U{z,z), defined 
on C or C*. The gradient of any real valued function of z and z is then given by : 

W(r)^2%l) (1) 

where^=^=i(^+z^). 

The equation of motion for our system takes the form: 

■■ r.dUiz,z) 

z + 2 = 2 

oz 

where the dot represents the time derivative. 

For two vectors A and B in the {0,u^,u^) plane, A x B is represented in complex form by the real quantity 
Im (AB) , A and B be being the complex affixes of A and B respectively. 

The angular momentum (which necesseraly conserves its direction orthogonal to the plane) L (t) = ^ (t) x ^ (<) = 
L{t) k , admits then the following correspondent : 

L ^Im{z{t)z{t)) ^ (zz ~zz^ (3) 
In the case of a central potential U{z,T) — U{\z\) — U{r), we have: 

VU{r) ^ -U'ir) (4) 
r 

where f'{x) = ^^^^ is the usual derivative of a function / of one variable x. 
Eq.([21) becomes simply: 

'z+-U'(r)=0 (5) 
r 

Using this result, the angular momentum's conservation is immediate (see Eq.Q): 

L= -^(jz-Jz^ =0 (6) 

2. Radial equation of motion 
If we use a polar representation for z, z — re'^, with 9 e [0, 2tt[, we have: 

(r^) = 2rr — zz + zJ (7) 

which gives: 
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and (since e^*^ = =) ' 



Then, we can write: 



and: 



z ^ re'" + lOre'" = [ r + i-\ - (10) 

r r 



^=(;-^)^ (11) 



Then using Eq.©, we deduce an equation for the the radial motion: 

L2 



which is readily integrated as : 



« =\E-^-U{r)j (13) 
i5 is a constant of integration, which can be identified with the energy of the system, since from Eq. PU)) we have: 

1|.|2 1/f . L\ f r . L \ _ 1//.n2 L 



Z\ hi^ Z l—^\Z 



If we note T{r) ^ E - U{r), Eq.([12l) takes the form: 



{ry + -\=E-Uir) (14) 



where 



fir) = '^T{r) 1 = ^{E F^r)), (16) 
^l(^) = U{r) + ^ being the radial effective potential. 

From Ea. (|15p and Eg. pop . we see that the radial oscillation motion presents, for bounded orbits, two turning point 
Tm and tm, corresponding to the pericenters et apocenters (apsidal positions) of the orbital motion. They are roots 
of the equation: 

f{r)^O^VLir)^E (17) 

In the sequel, we will suppose that Vm and vm, are simple roots, that is V[{rm) < and V^(7'm) < 0. 
At the apsidal positions, r,„ = vm — and: 

~ ^^j^j ~ *~2'^i' j — rn, M (18) 
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The i factor in the right hand side of Eg. ljTS)) imphes that, for apsidal positions, the instantaneous speed vector is 
orthogonal to the position vector. 

Note that if z(t) and z{t) are continuous and differentiable for every values of t, the differentiability domain of r{t) 
is stricly limited to the real interval ]rm,rj\/[. 

To extract from Eq. (fT5)) the radial celerity r{t) as a real- valued function, the square root has to be constructed 
with a pre-defined sign. To take account of the turn back in the radial motion, this sign changes when r reaches the 
extremal values and rjv/, which corresponds to change the square root determination. 

We obtain: 



L 



■^[t)=± — JWW) (19) 
r{t) 



or: 



dr{t)^±-^J{r)dt (20) 
r 

The -I- sign is chosen when r moves from to tm (increasing phase of the radial oscillation, dr{t) > 0) and the 
- sign during the reverse motion (decreasing phase of the radial oscillation, dr{t) < 0). To be more precise, we will 
index each phase of the motion by an integer number fc e N, the even value k — 2n corresponding to the decreasing 
phase of the n-th period (except for n = 0) and the odd value k — 2n + 1 corresponding to the increasing phase of 
the same period. With this convention, we then have: 

drit) = {-l)'^+'-^^/mt))dt (21) 

and: 

L 



m^i-rf+'^^^mt)) (22) 



Eq. ([2T|) gives implicitely the solution of the radial equation of motion Eq. p2 



3. Orbital equation and r-parametrization of the motion 



In the fc-phase of the motion ,in every point except the apsidal ones, r{t) is a monotonic function (see Eg. ([22]) ) and 
we can choose to parametrize the motion by r rather than by t. 
For instance: 

'^^ ^ ' -1)'^+^^= (23) 



drit) r r/7W 

where we have used Ea. (p2)) and Eq.®. 

Choosing a reference time to during the same phase, we can write: 



>{r{t)) = e{ro) + {-ir+' (24) 



with r(io) — fo and the orbital equation for this phase writes: 



e{r)^e{ro) + {~lf+' (25) 

We have to be more careful if we want to use r as a global parameter for the motion. Indeed, since r{t) is a periodic 
bmction of t, t{r) and every monotonic function of as 9 (see Eq.®), will be multivalued. Therefore, at every k 
phase of the motion will be attached a different branch of the function. 
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For simplicity, we choose the initial condition: 

" r{t = 0) = r„. 

e{t = o) = o 

Then, the different determinations of 6 (r) corresponding to each phase of the motion are given by: 



(26) 



where: 



In the first phase of the motion : 9 {r{t)) — g{r{t)) (27) 

dp 

In the second phase of the motion : 9 {r(t)) = $ — / — j= = 2$ — g{r{t)) 
In the third phase of the motion : 9 {r{t)) = 2$ + g{r{t)) 

(28) 



/"■f*) do 

g{r{t)) = I (29) 



$ = g(fAf ) = ll" r'^^'fir') ^^^^^ apsidal angle. 
More generally, we will write: 

In the fc-th phase of the motion : 9k (r(t)) = 2n$ + (-l)'=+\g (r(i)) (30) 
where n is the integral part of ^. 

From Eq. ([22|) and Eq. ([23|l . we deduce the following expression for the instantaneous velocity in the fc-th phase: 

m = {i-lf^'^/mt))+i) = A ^/fm)+^) e^'"^^^'^^ (31) 

9k {r) being given by Eq.pO|). 

Note that, since z = \z where A has a nonzero imaginary part, Eq. (j3ip ensures the non-colinearity of and ~r*. 
In the sequel, we will have to use the vector x L which has the following complex correspondent : 

y (i) X L ^ zLz{t) = - (l + (-1)' ^^mt))) (32) 
We see immediately that i x and are always linearly independent, except at the apsidal positions. 

III. PERES APPROACH TO THE FBRS VECTOR 

In the special case of the Kepler system U{r) ~ — r^^ , the system admits the Laplace-Runge-Lenz vector, Alrl 
— ^ X L — as a supplementary conserved quantity^. In the complex formulation used above, this gives : 

Alrl Alrl = iLz - -z (33) 

r 

Following Peresi^, for the general central potential U{r), we will then look for a supplementary conserved vector of 
the form: 



A^r X L [ jj^air) ] + b{r)r (34) 
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whose associated complex correspondent is: 



y.2 



A = + KO^, (35) 

a(r) and 6(r) being two real valued functions that we will have to determine (the factor in the first term has 
been introduced for future convenience). 

Using Eg. ([32|) . we can rewrite A in the k-ih phase, as: 

Ak = {b{r) - a{r) + {-lf+^ i^/J{^a(r)) z (36) 

4- Differential system for the coefficients of the FBRS vector 
On every time interval / on which A is constant we must have: 

^ = (37) 

that is, with Eq. ipTjl : 



= I {h{r) - a{r) + (-1)'=+^ ^VWHr)) z + {h{r) - a{r) + *v^«(r)) z (38) 

r (^b'ir) a'{r) + {-lf+' ^ i^^)a'{r) + ^^)=) ) ^ 

+ {h{r) ~ a{r) + {-l)'"-' iy/W)a{r)) ((-1)'+' ^/W) + = (39) 



Using Eq. iP^ . this becomes: 

i-lf^'^) {h'{r) - a'[r) + + , i^^)a'{r) + (^J^^ + l^^^ a(r) + b{r)^ = (40) 

Consequently, in the fc-phase of the motion, on the interval of radial values ]rm,rM[, we must have: 

r a'(^)_6'(r) + Ml)_Ki) ^0 

I rf{r)a'{r) + a{r) (r^^ + f{r) - l) + b{r) = ^^^^ 
Extracting a'{r) from the second equation above and substituting in the first one, we obtain, in matrix form: 

^a{r)\' _ -1 ( M + f(r)-l 1 \ fa{ry 



h{r)J rj{r) \ r^li - /(r) - 1 /(r) + 1 / \h{r) 



(42) 



5. Exact solution 

Peres^, starting from a first order differential system equivalent to Eq. ([4T|) transforms it in a second order differential 
equation for the coefficient a(r). Nevertheless, the structure of this last equation is, at first sight, rather complicated 
and the author restricts its analysis to a characterization of the singularities. 

Things are much more transparent if we choose to work with a slightly different unknown function. 
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Indeed, it is readily seen that the first equation in Eq. pT]) can be rewritten: 

m'(r) + 2a{r) = (r6(r))' (43) 
Introducing an auxihary function u{r) defined by: 

a{r)=u'{r), (44) 

the above equation Eq. (|4T]) becomes: 



ru"{r) + 2u'{r) = {ru{r))" = {rbir))' (45) 

Then: 

h{r) = - (ru(r))' + &o = u'(r) + -u(r) (46) 
r r 

where the integration constant has been chosen equal to 0. 

We can report the expressions Eq. (j44|) and Eq. (j46|l of a{r) and h{r) in terms of u{r) in the second equation of 
Eq.dH]). This gives: 

(u'{r) + \u{r)^ = -j^ (^u'ir) (^r^ /(r) - + {u' (r) + i«(r)) (/(r) + 1)) (47) 



or 



1 f'jr) \ 1 
r + 2/(r)r^'^ + r2/(r) 



+ ( r + TTFH^ ] ^'(r) + ^^^/(r) = (48) 



In a more compact form, we obtain the following second order differential equation for u(r): 



2 



or (see Eq.^): 



""(0 - ^«'(^) + (g'iryf u{r) = (50) 
The coefficients in Eg . ([501 depending only on g (r), if we define: 

u{r)^v{g{r)) (51) 

Eg. ([50]) becomes then: 

[g' {r)f v^'Hg (r)) + g" (r) v^^\g (r)) - ^ (g' (r) v^^\g (r))) + {g' {r)f v{g (r)) = (52) 
that is, simply: 



v^^\g {r))^v{g (r)) = 



(53) 
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where v^"\x) = Qk- 

The resolution of Eg. ([55)) is immediate and gives: 

u(r) = Q!Cos.g (r) + /3sing (r) , Q!,/3eR (54) 
Reporting this resuh in Ea. (|^^ and Eq. (|^5)l . wc finally have: 

f °(^) ^ A, , (/? cos g(^) -Q sin g(r)) 

" ^/('■) ^ , a,/? e M (55) 



that is: 



6(r) = a(r) + i (a cos 5 (r) + /3 sin g (r) ) 



a(r) — — 3= cos (g (r) + (h) 

" ^\ . ' ' , 7eK, 0e [0,27r[ (56) 



1^ 6(r) = a{r) + ^ sin {g (r) + 1 

6. Discontinuous behavior of A 

The preceding result Eq. (j56p determines completely the form of A on ]r„i, tmI during a given phase of the motion. 
Indeed, inserting Eq. ([56l) in Ea. ((36)) . we obtain: 



Ak = J (sin {g [r) + 0) + {-if+H cos (5 (r) + 0)) z{r) = (-l)fc+ii7e(-i)'*(9W+«e*«^M (57) 
Using Eg. ipO)) . this becomes: 

A = (-l)'^+ii7e(~^)'°*(3(''^+'^^e*(^"*+("^^''^'f(''^) = (-l)'=+i'yie("^)'°"^+2™* (58) 

where n is the integer part of ^. 
Under a more detailed form, we have: 



-42,1 = -7e*f +"^+2™* 

A2n+l = 76*5^*^+2™* 



(59) 



The integration constants 7 and 4> for and determine then respectively the modulus and argument of the A vector 
afhx. The freedom in the choice of the parameters 7 and cf) induces that A can be identified with any vector of 
the plane. We have the same type of result for the specific Kepler problem. Indeed, in this case, we can add to 
the Laplace-Runge-Lenz vector Ak a second conserved vector, S k = L x Ak, the so-called Hamilton vectori^iii, 
whose constancy is a direct consequence of these ones of L and Ak- Any linear combination of this two orthogonal 
vectors being conserved, we can build a conserved vector corresponding to any vector of the plane. S k and Ak are 
particularly interesting choices because they give the directions of the minor and major axes of the elliptical trajectory. 

If we want that the functional forms of a(r) and 6(r) being globally defined during all the motion, we have to keep 
the same values for 7 and (j) in every phase. With the choice (p — ^, Eq. ((59| becomes: 

A2n = A2n+1 = 76^"^* (60) 

A is then constant when we pass from a even phase to the following odd one, that is when we cross the apocenter. 
But after a complete period of oscillation, when we reach the pericenter again, we pass from k = 2n to k ~ 2n + I. 
At this moment the value of .42ti+i change toy^2n+i with: 

A2n+2 = A„+ie''* (61) 

This corresponds to a 2$ rotation of the associated vector. We recover here the discontinuity jumps (observed by 
Serebrennikov, Shabad, Buch and Denma n^^i^^ and Peres^^ and first studied in a detailed by Holas and March^^) 
which make FBRS vector only a piecewise conserved quantity: the perihelion vector A presents discontinuities at 
each pericenter (with the choice made here and at each apocenter if we choose the initial condition r{t = 0) — tm) 
corresponding to a rotation of two times the apsidal angle. 
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7. Harmonic oscillator, Kepler problem and Bertrand's theorem 

For $ = TT, which is the case of the Kepler problem U{r) = — r^^ for every values of the characteristic parameters 
of the motion L and E, A is a true vector conserved quantity and is identical to the usual Laplace-Runge-Lenz vector 
Ak- Indeed, in this case, with E < 0, introducing Clairaut's variable* u — ^, we have (see Eq. ([29| and(see Eq. (fl6|) ): 



= 7k 



L 

\/2 Jti,„ = - 



dv 



which givea^ for 1 - 2L^ \E\ = > OS: 



argsm 



\E\ 



2 



(62) 



(63) 



Since is the smallest root of /(r) = (^—\E\ + \ — — 0, that is, Um = 7- is the greatest root of 



\E\+u — k^v? = 0, we obtain: 



1 + e 



(64) 



and: 



Then (see Eq.(l56l) with = f ) 



g{r) =^ argsin 



1-^ 



TT 

2 



(65) 



a{r) 



:\/i-(^r 



-7: 



b{r) = a{r) + 2 sin {g (r) + f ) = -7^ " ? ^ = 
Taking 7 = — e, we recover the coefficients of the usual Laplace-Runge-Lenz vector (see Eq. 



-7- 



(66) 



a(r) 
6(r) 



(67) 



In the isotropic harmonic oscillator problem (Hooke's problem) $ = ^, for every L and E. The global direction is 
then conserved but the sense of A is alterned at each pericenter crossing : A ^ — li^. The generalized Hamilton 
vector S = L y. Ais subject to the same phenomenon : S ~> — S . for the Fradkin tensor T = ®^ + ®^ 



2\A\ 



■S ® S ^ we see immediately that it's a global invariant 



11 — — 

, recalling that it can be written as T = ^ A ® A 
of the motion. 

As established by Bertran d^^i^^ more than one century ago, Kepler and Hooke potentials are the only central 
potentials for which the apsidal angle is commensurable with tt for every values of the initial parameters of the 
motion. For all the other central potentials, this condition, which is necessary for the closure of the orbit, is obtained 
only for specific values of E and L. In these cases, as established by Holas and MarchSi^ it is still possible for an orbit 
of multiplicity n to build global geometrical invariants in form of n-arm stars by using the n distinct FBRS vectors 
associated to the system. Nevertheless, as we have seen before, the existence of a general (that is for every initial 
condition) true invariant vector or tensor is a specificity of Kepler and Hooke problems respectively. 
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